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Abstract 



An updated search for heavy narrow resonances decaying to muon or electron pairs 
using the CMS detector is presented. Data samples from pp collisions at a/s = 7 TeV 
and 8 TeV at the LHC, with integrated luminosities of up to 5.3 and 4.1 fb _1 , respec- 
tively are combined. No evidence for a heavy narrow resonance is observed. The 
analysis of the combined data sets excludes, at 95% confidence level, a Sequential 
Standard Model Zg SM resonance lighter than 2590 GeV, a superstring-inspired ZL 
lighter than 2260 GeV, and Kaluza-Klein gravitons lighter than 2390 (2030) GeV, as- 
suming that the coupling parameter k/ Mpi is 0.10 (0.05). These are the most stringent 
limits to date. 
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1 Introduction 

A number of scenarios for physics beyond the standard model predict the existence of heavy 
narrow resonances that decay to lepton pairs. In this Letter, we report on a search for reso- 
nances with the Compact Muon Solenoid (CMS) detector at the CERN Large Hadron Collider 
(LHC) HI. We consider the following three benchmark scenarios: the Sequential Standard 
Model Zg SM with standard-model-like couplings [2J, the Z'^ predicted by grand unified theo- 
ries 0, and Kaluza-Klein graviton excitations in the Randall-Sundrum (RS) model of extra 
dimensions |3H5|. The RS model has two free parameters. One parameter is the mass of the 
first graviton excitation, and the other is the coupling k/ Mpi, where k is the curvature of the 
extra dimension and Mpi is the reduced Planck scale. 

Previous searches for narrow Z' — » £ + £~(£ = }i,e) resonances have been reported by the 
CMS [6] and ATLAS [7] Collaborations, each based on integrated luminosities of 5 fb _1 at y/s = 
7 TeV. The CDF and DO experiments have published results based on integrated luminosities 
exceeding 5fb _1 of pp collisions at y/s = 1.96 TeV H8HT3H. The best previous direct lower limits 
on the Zg SM and Z'^ masses are 2330 GeV and 2000 GeV [6[, respectively. The best previous 
direct limits on RS graviton (Gkk) production are 2160 GeV for k/M?\ = 0.1 [7] and 1810 GeV 
for k/ Mpi = 0.05 [6|. Indirect constraints IfMHTTI are less stringent. 

We use data samples from pp collisions at center-of-mass energy y/s = 8 TeV, corresponding to 
an integrated luminosity of 3.6 ± 0.2 fb _1 for the dielectron channel. The dimuon channel does 
not use information from the calorimeters and incorporates additional data from running peri- 
ods when the calorimeters were not fully operational. This increases the integrated luminosity 
of the sample to 4.1 ± 0.2 fb -1 . We combine the analysis of these data with previous results 
from the analysis based on an integrated luminosity of 5fb _1 at y/s = 7 TeV J5|. The recon- 
struction, selection criteria, efficiencies, and systematics for the two data sets are very similar. 
The results are applicable to any model with a narrow resonance that has equal dimuon and 
dielectron branching fractions. A resonance is considered narrow if the experimental width is 
dominated by the detector resolution. 

We perform a shape-based analysis of the dilepton mass spectra searching for a peak on a 
smoothly falling distribution with the overall background normalization determined by an 
unbinned maximum likelihood fit. The data are consistent with expectations from the standard 
model. We report limits on the ratio (R^) of the production cross sections times branching 
fractions of a heavy narrow resonance to that of the Z boson, at the 95% confidence level (CL). 
Many experimental and theoretical uncertainties cancel in this ratio. We further translate these 
limits into lower limits on the masses of new heavy narrow resonances, using next-to-next- 
leading-order (NNLO) cross section calculations [18] for the Z boson production. 

2 The CMS detector 

A detailed description of the CMS detector can be found in Ref. Ifl9l . We briefly discuss the sys- 
tems most relevant to this analysis. The central feature of the CMS detector is an all-silicon in- 
ner tracker system, composed of silicon pixel and strip detectors. The tracker is surrounded by 
a lead-tungstate crystal electromagnetic calorimeter (ECAL) and a brass-scintillator hadron cal- 
orimeter (HCAL). The finely segmented ECAL consists of nearly 76000 lead-tungstate crystals 
which provide coverage up to pseudorapidity \fj\ = 3.0. It is divided in the barrel (\rj\ < 1.479) 
and endcap (1.479 < \t]\ < 3.0) detectors. We define pseudorapidity as rj = — In [tan (6/2)]. 
Here, 6 is the polar angle with respect to the direction of the counterclockwise proton beam. 
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3 Event selection and object reconstruction 



We use <p for the azimuthal angle of a track's momentum at the point of closest approach to the 
beamline. The tracker and calorimeter systems reside within a 6 m diameter superconducting 
solenoid, which produces a 3.8 T axial magnetic field. Muons are detected by gas-ionization 
chambers embedded in the steel flux-return yoke. 

The CMS experiment utilizes a two-level trigger system. The first level of the trigger (LI) selects 
events of interest using custom hardware processors |20|. It uses information from the muon 
and calorimeter systems to reduce the readout rate from the 20 MHz bunch crossing rate to 
a maximum rate of 100 kHz. The software based high-level trigger (HLT) further reduces the 
recorded event rate to a few hundred Hz by adding information from the inner tracker and 
analyzing event information in greater detail [21 j. 

3 Event selection and object reconstruction 

The event selection closely mirrors the one used for the \/s = 7 TeV analysis. We briefly re- 
view the procedure here. Dimuon events are triggered by requiring at least one muon to be 
reconstructed by the HLT and to have transverse momentum pj > 40GeV and \rj\ < 2.1. Di- 
electron events are accepted by a double-electron trigger requiring two clusters in the ECAL, 
each with transverse energy Ej > 33 GeV. The trigger allows only small deposits of energy 
in the HCAL to be associated with the ECAL clusters. HLT clusters are required to be loosely 
matched to the trajectories of tracks having hits in the pixel detector. The lepton trigger effi- 
ciencies are measured using a "tag-and-probe" technique at the Z resonance j6j|22H231, up to 
transverse momenta of roughly 500 GeV for muons and 100 GeV for electrons. For higher trans- 
verse momenta, the electron efficiency is measured using a simple trigger which requires only 
an ECAL cluster with Ej > 300 GeV to directly monitor the trigger efficiency for selected high 
mass events. In order to have a consistent trigger between the low-mass control region and the 
high-mass signal region, the simple ECAL trigger is used only to validate the efficiency of the 
primary trigger. The muon trigger efficiency is 97% for events with both muons within the trig- 
ger acceptance, across the entire range of dimuon invariant masses of interest. The efficiency 
of the electron trigger for dielectron candidates passing the analysis selection requirements in- 
creases from 80% at electron Ej = 35 GeV to a 99% plateau at Ej > 37 GeV. The efficiency 
threshold curve of the trigger is measured using data collected by a lower threshold trigger 
that is applied to approximately every 5 th event passing the LI part of the trigger. Because 
the threshold behavior is well-determined, the offline Ej selection cut can be placed at 35 GeV, 
which improves the normalization to the Z peak. Both muon and electron trigger efficiencies 
are within 1-2% of those found in Ref. |6j. Standard CMS algorithms I|6l l2"3"ll2"4| are used to re- 
construct and select muon and electron candidates. Muon candidates are formed by matching 
tracks in the silicon tracker to tracks in the muon systems. Muon tracks are required to have 
hits in nine or more layers of the tracker and include at least one hit from each of the pixel and 
muon systems. 

Muon candidates are required to be isolated in a cone about the muon direction of AR < 0.3 
in the tracker and to have pj > 45 GeV and \tj\ < 2.4. The quantity AR is defined as AR = 
a/ (A^) 2 + (A(p) 2 , where A(p is in radians. The combined fit of the muon trajectory through 
the tracker and muon systems provides a reliable measurement of muon momenta extending 
to order 1 TeV p3l l25| . Electron candidates are formed by matching ECAL clusters to recon- 
structed tracks, and are required to have \rj\ < 1.442 and 1.560 < \t]\ < 2.5 for the ECAL barrel 
and ECAL endcap regions, respectively. As in Ref. [6J, electrons are additionally required to 
have little associated activity in the HCAL, have a shower shape consistent with that of an 
electromagnetic object, and be isolated in a cone about the electron direction of AR < 0.3 in 
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the calorimeter and tracker. To account for the contamination which varies with the number 
of additional proton-proton interactions per event, the calorimeter isolation of the electron is 
corrected for the average energy density in the event [26 j. 

The ECAL is capable of measuring energies in a single crystal of up to approximately 1.7 TeV 
in the barrel region and 2.8 TeV in the endcap region. Saturation effects become an issue for 
resonance masses around 4 TeV, which are beyond the reach of the current analysis, within 
the framework of the models considered here. Occasionally, anomalously large signals are 
observed in the barrel ECAL due to the direct deposition of energy into the avalanche photo- 
diodes (APDs) by particles transiting the detector p7[. Since the APDs are normally used to 
detect the scintillation light produced by the crystals, the equivalent energies of these signals 
can reach into the TeV range. The deposits are generally in a single channel and have different 
pulse timings from scintillation signals. They are rejected by cutting on the timing of the pulse 
and the amount of energy recorded in neighboring crystals. 

More than 99% of the anomalous clusters are rejected with a negligible loss of real electrons [27J. 
When combined with the requirement of a compatible track in the silicon tracker, this require- 
ment reduces background from anomalous clusters to a negligible level. However, the energy 
deposit pattern of all selected high energy electrons is scrutinized even further. There is no 
evidence that any of the remaining events is due to these anomalous signals. 

Previously established techniques |22l [231 are used to measure event reconstruction and se- 
lection efficiencies. The efficiencies for reconstructing and selecting leptons are roughly 90%, 
for both muons and electrons with pj > 100 GeV. The efficiencies include the isolation re- 
quirements and are measured with total uncertainties of a few percent. Part of the systematic 
uncertainty cancels in the ratio of high-mass dilepton cross section to the Z boson cross section. 
Studies with Monte Carlo simulated samples predict both lepton reconstruction and selection 
efficiencies to be constant within 1% for transverse momenta pj > 100 GeV. 

Each dilepton candidate is required to have two isolated leptons of the same flavor that pass 
the identification criteria described above. When multiple dilepton candidates are present, the 
most energetic pair in the bunch crossing is selected. For dimuon events, one muon must have 
\t] | < 2.1 to satisfy the trigger requirements. For dielectron events, at least one electron must 
have \tj \ < 1.442. This removes events in which both electrons are in the endcap, a topology 
where little signal is expected but which has a significant background arising from misidenti- 
fied jets. Muons are required to have opposite charge, since a charge mis-assignment implies 
a large mismeasurement of momentum. The energy estimate for electrons is dominated by 
electromagnetic calorimeter information and is not sensitive the momentum mis-measurement 
indicated by a charge mis-assignment. The charge requirement would also result in a few per- 
cent efficiency loss in a region with little background and would degrade the sensitivity of the 
analysis. Therefore, we do not impose a charge requirement for dielectron candidates at high 
mass. Muon candidates are additionally required to originate from the same vertex. The x 2 
per degree of freedom for the fit to a common vertex is required to be less than 10. The tracker- 
measured transverse impact parameter with respect to the beam spot must be less than 2 mm 
for each muon. The mass resolution of a dielectron candidate is predicted by Monte Carlo sim- 
ulation to be approximately 1.8% for masses above 800 GeV. The dimuon mass resolution is 5% 
(9%) at 1 TeV (2 TeV) and increases linearly with dimuon mass. 

The opening angle of the muon pair is required to be less than n — 0.02 radians. This re- 
quirement greatly reduces the cosmic ray background associated with muons traversing the 
detector. 
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5 Results 



4 Backgrounds 

The dilepton background in pp collisions at \/s = 8 TeV is very similar to that found in y/s = 
7 TeV collisions |6l even though there were significantly more interactions per bunch crossing 
at the higher energy. The effect of this "event pileup" is included in our simulations of back- 
ground processes and our data-driven estimates of the background from misidentified jets. The 
dominant standard model background is due to Drell-Yan production. The shape of this contri- 
bution is determined from Monte Carlo simulation using the pythia v6.4 [28] event generator. 
The background contribution is normalized to the event count at the Z peak by counting same- 
flavor dilepton candidates within the mass window 60 < < 120 GeV. The next largest 
background contribution is due to other standard model processes that produce isolated dilep- 
tons. We consider the lepton flavor symmetric processes of tt, tW, Z — > t + t~, and diboson 
(WW, WZ, and ZZ) production when estimating this background component. The absolute 
normalization and shape for these backgrounds is taken directly from Monte Carlo simulation 
generated using MadGraph 5 [29], POWHEG Il30ti33| , and PYTHIA. We validate this back- 
ground prediction by comparing the e}i dilepton mass spectra for data and simulation. 

Track-based lepton isolation strongly suppresses backgrounds from jets misreconstructed as 
leptons. This background is almost negligible for the dimuon channel but is a significant por- 
tion of the non Drell-Yan background in the dielectron channel. Since misidentification of jets 
as leptons is more likely to occur for electrons than for muons, electrons have additional iso- 
lation and jet discrimination requirements. In the dielectron channel, the main contributing 
processes apart from Drell-Yan are dijet, W+jet, and photon+jet production (where the photon 
is misidentified as an electron). The probability that a jet is misidentified as an electron is mea- 
sured in bins of Ej and t], using a jet-dominated sample. This probability is then used to weight 
events in which one electron satisfies all selection criteria and the other is a candidate for being 
a misidentified jet, to obtain the jet background prediction. The dimuon resonance search is 
susceptible to backgrounds from cosmic ray muons. The expected cosmic ray background for 
dimuons with m fl{l > 200 GeV is determined from two complementary samples. For events 
in the first sample, the requirement on the dimuon opening angle is removed. In the second 
sample, the impact parameter requirement on the muon tracks is not applied. From the popu- 
lations of these two samples, the remaining cosmic ray background contamination is estimated 
to be less than 0.2 events. 

5 Results 

The dilepton mass distributions for events passing all the selection criteria are shown in Fig.JT] 
The "jets" distribution illustrates the contribution of events in which at least one jet is misre- 
constructed as a lepton. This distribution is derived from data while all other components are 
derived from simulation. The relative fractions of the different background components are 
fixed to the ratios of their theoretical cross sections. The total simulated background is normal- 
ized to data at the Z peak (60 < mm < 120 GeV). The expected yields in the control region 
(120 < nift > < 200 GeV) and in the search region (ma > 200 GeV) are compared with observed 
yields in Table [I] The observed mass spectra and event counts agree with standard model 
predictions both in shape and normalization. 

We set a 95% CL limit on the ratio R a of the product of the cross section and branching fraction 
for each Z' boson to that of the standard model Z boson. The cross section of the 21 boson 
is calculated in a window of ±40% about the on-shell mass of the resonance, while for the Z 
boson it is calculated in the peak window defined above. We follow the Bayesian procedure of 
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m(n+n") [GeV] m(ee) [GeV] 



Figure 1: The invariant mass spectrum of (left) and ee (right) events for the = 8 TeV 

data set. The points with error bars represent the data. The solid histograms represent the 
standard model predicted background contributions. 

Table 1: The dilepton event count in the control region 120 < > < 200 GeV and in the search 
region mu > 200 GeV for the \/s = 8 TeV data set. The total background is the sum of the 
events for the standard model processes listed. Uncertainties represent a quadrature sum of 
statistical and systematic uncertainties. 



Source 




Number of events 






Dimuon sample 


Dielectron sample 




(120-200) GeV 


>200GeV 


(120-200) GeV 


>200GeV 


Data 


1381 


3503 


12030 


2904 


Total background 


13010 ± 590 


3630 ± 160 


12240 ± 590 


2970 ± 260 


Z/7* 


11700 ±570 


2920 ± 140 


10660 ± 530 


2200 ± 220 


tt + others 


1280 ± 150 


698 ± 78 


1220 ± 180 


560 ± 80 


Jets 


26 ±3 


10 ±1 


360 ± 180 


210 ±110 



Ref. |6|, which is based on an unbinned extended maximum likelihood analysis. We calculate 
the limits using the 8 TeV data alone, as well as from a combination of the 8 TeV and 7 TeV data 
sets. Mass-dependent ratios of parton distribution functions (PDF) at y/s = 7 TeV and 8 TeV 
are used as an additional input to derive limits on R a at 8 TeV, R^STeV, that combine both data 
sets. The CTEQ6.1 LO PDF set |34| was used to calculate these ratios, and the result was cross- 
checked with the MSTW2008 PDF set (351 . The CTEQ and MSTW calculations agreed well and 
the uncertainty in this ratio does not significantly contribute to the final result. The most signif- 
icant uncertainty in the limit computation is associated with our understanding of the selection 
efficiency and detector acceptance ratio for Z' bosons relative to the Z, denoted R e . The uncer- 
tainty in the total lepton selection efficiency at high mass dominates the R e uncertainty. The 
lepton selection efficiencies are measured in data up to pj ~ 500 GeV, but above 100 GeV the 
uncertainties in these measurements become large. This leads to a total uncertainty in R £ of 3% 
for the dimuon channel and 8% for the dielectron channel after including PDF uncertainties in 
the acceptance. The effects of misalignment, higher order corrections to the background shape, 
and the uncertainty in backgrounds due to jets misidentified as leptons have only negligible 
impact on the limits. The upper limits on the ratio R a for spin-1 and spin-2 particles obtained 
from the dilepton combined mass spectra are shown in Fig. [2] Table [2] shows the limits on R a 
converted into mass limits on specific models. The resonance is assumed to be narrow, mean- 
ing that the detector resolution dominates the width of the peak. The with a relative width 
of 0.6% is therefore considered narrow. A wider resonance, such as the Zcc M , which has a width 
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6 Summary 



of 3%, will have more background under the peak. Consequently, we would set weaker limits 
on its production cross sections. The two cases provide similar results when there is very little 
background after all selection criteria have been imposed. This occurs around 1.4 TeV. For a res- 
onance below 1.4 TeV not to have been discovered, it must have a small coupling and therefore 
be narrow. For the spin-2 case an additional requirement is that the ratio of gg to qq production 
of the resonance must be the same as the ratio for an RS graviton. The combination of the 7 and 
8 TeV data sets relies on this assumption, as gg and qq cross sections scale differently with \fs. 
For the spin-1 case, no gg coupling is considered. The Z' and RS Graviton cross sections are 
calculated using the PYTHIA event generator with the CTEQ6.1 PDF set. The LO cross sections 
are corrected for next-to leading (NLO) or NNLO QCD contributions using the same k-factors 
as Ref. 13. A mass dependent NNLO k-factor calculated with zwprodp ||36ti38| is used for the 
Z' models. A flat NLO k-factor of 1.6 is applied to the RS graviton cross sections |39J. 
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Figure 2: Upper limits on the ratio R a of the production cross section times branching fraction 
into lepton pairs to the same quantity for Z bosons, as a function of resonance mass M for spin- 
1 (top) and spin-2 (bottom) boson production. The left plots are for the 8 TeV data set while 
the right plots are for the combination of the 7 and 8 TeV data sets. For the spin-2 case, the 
7 and 8 TeV data set combination is only valid for models that have the same fraction of qq 
to gg coupling as an RS graviton. For the spin-1 case no gg coupling is considered. Shaded 
bands identified in the legend correspond to the 68% and 95% quantiles for the expected limits, 
respectively. 



6 Summary 

The CMS Collaboration has searched for heavy narrow resonances in dimuon and dielectron 
invariant mass spectra. The search combined data samples from pp collisions at a/s = 7 TeV HQ 
and 8 TeV. The \/s = 8 TeV data sets have integrated luminosities of 4.1 fb -1 (3.6 fb -1 ) for 
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Table 2: Mass lower limits at the 95% CL on specific models obtained using dilepton data at 
y/s = 7 and 8TeV separately and combined. The 7TeV results are taken from Ref. [6|. 



Model 


Mass Limits (GeV) 


7TeV 


8TeV 


7+8 TeV 


7' 




2330 


2440 


2590 


7' 




2000 


2110 


2270 


G KK {k/M n 


= 0.1) 


2140 


2260 


2390 


Gkk (k/M F1 


= 0.05) 


1810 


1900 


2030 



the dimuon (dielectron) channel. The \/s = 7 TeV data sets have integrated luminosities of 
5.3 fb _1 (5.0 fb _1 ) for the dimuon (dielectron) channel, and have been previously published 0. 
The measured dilepton mass spectra are consistent with predictions from the standard model. 
Upper limits on the cross section times branching fraction for the production of new heavy 
narrow resonances relative to Z boson production are presented. The findings exclude, at 95% 
CL, a Zg SM with standard model-like couplings below 2590 GeV and the superstring-inspired 
below 2260 GeV. An RS graviton with k/Mn of 0.1 (0.05) is excluded below 2390 (2030) GeV. 
These are the most restrictive limits to date for the classes of models considered. 
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